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We calculate the contribution of the two photon production process into e^ e~ spectra, and com- 
pare the results with experimental data from the PHENIX detector at RHIC. We study the contribu- 
tion given by "semi-coherent" kinematics, in which one photon is relatively hard and is incoherently 
emitted by participating protons, while another can be soft enough to be in a coherent domain. 
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I. INTRODUCTION 

Thirty years ago one of us [D [2] had suggested to use 
dileptons and photons as "penetrating probes" for dense 
hadronic matter created in ultrarelativisitc heavy ion col- 
lisions, which - unlike hadrons - are observable from all 
stages of the collisions and thus can tell us what the ini- 
tial hottest temperature reached can be. It is a very 
challenging task for experiments, as one has to remove 
hadronic backgrounds orders of magnitude larger than 
the photon or dilepton signal. And yet, over the years 
there were successful measurements, both at CERN SPS 
(muon pairs by NA50/NA60, electron pairs in CERES, 
photons in WA98) and RHIC (photons, muons and elec- 
trons in PHENIX, electrons in STAR). We will not go 
into details of these works, just make few general com- 
ments. 

Already the above mentioned papers from 1970's have 
singled out the so called intermediate mass dileptons 
(IMD's), with the mass 1-3 GeV or between </> and J Ii\j 
resonances, as the window for observing the thermal 
QGP radiation. More detailed predictions have been 
made in Ref. |4J, where it has also been predicted that 
most of those pairs observed are not from charm decays, 
as was widely believed at the time. Only with success- 
ful completion of the NA60 experiment, with its sophis- 
ticated charm tracking, this collaboration had recently 
confirmed that they do indeed observe thermal radiation 
from QGP 5 and not just charm decays. For summary 
of other NA60 results see e.g.0: those include dileptons 
with small masses which come from resonances p mesons 
decaying in hadronic and near-Tc region. Although still 
far from being perfect, the existing theory provides a rea- 
sonable overall description of the NA60, see e.g. [21 [3]. 
Important recent observation of thermal photon radia- 
tion from hadronic gas and QGP has been also made by 
PHENIX collaboration ^6], which is also in overall fair 
agreement with the current theory and the hydrodynam- 
ical picture of the collision. 

And yet, some aspects of the experimental data at 
RHIC remain puzzling. Dilepton results from PHENIX 
show production rate of small mass M ~ 500 MeV dilep- 
tons few times above theory predictions. Another puzzle 
is the presence of the so called "cold" component in the 
dilepton spectrum for pt < 500 MeV, which is shown in 
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FIG. 1: (Color online) Acceptance corrected invariant e^e~ 
yield versus total transverse momentum of the dilepton pair, 
for pp collisions (left) and AuAu collisions (right), from 
PHENIX publication [TU]. The solid curves show the expec- 
tation from the sum of the so called hadronic cocktail contri- 
bution plus charm decays. 



FigjT] While the pp data (points in the l.h.s.) agree 
rather well with "hadronic cocktail" (curves), in AuAu 
data (r.h.s.) one finds systematic upward deviations of 
the data from from similar curves, at small pt (the left 
side of the AuAu plot). If fitted with exponential, the 
data have a slope Tg// ~ 100 MeV, which is about twice 
smaller than the typical slope of the main "hot" com- 
ponent. What is especially strange about it is that this 
slope seems to be the same for different dilepton mass 
bins, see the three lowest curves on the right hand side in. 
This is in contrast to the "hot" component, which shows 
Teff increasing with M, in good agreement with expec- 
tations based on hydrodynamical picture of expanding 
matter. It is a presence of such "cold" component which 
originally motivated us to have a look at some dilepton 
production mechanisms which are not included in the 
"standard" theory toolbox. 

Small mass component is another puzzle, it has un- 
usual centrality dependence. 

Coherent two-photon processes are a well-known 
source of small mass and small pt dileptons. Their basic 
theory had been developed already in 1930's, when Lan- 



dau and Lifshitz [TT] calculated their total cross section 
for nuclear collisions using Weizsacker- Williams (WW) 
approximation. There have been extensive studies of the 
so called ultraperipheral processes in RHIC environment, 
for experimental results from STAR collaboration see 
[l5] . As the name suggests, those processes take place 
at very large impact parameters b > 2R, at which no 
nuclear interactions take place. For electron pair pro- 
duction the characteristic b are related to the electron 
mass, and are thus very high. Theory development in- 
cluding all orders in Za has been worked out in the last 
decade. 

However the contribution of such processes at near- 
central collisions (when multiple hadronic production 
does happen) and for the kinematical range oi pt, M seen 
by PHENIX and NA60 has not to our knowledge been 
considered. This is what we are going to do in this work. 

Additional motivation for looking at the two-photon 
processes comes from the standard relations between on- 
shell and slightly virtual photons 7*, which are seen as 
small-mass dileptons. PHENIX has used such relations, 
relating dileptons with masses M > 100 MeV with real 
photons. However, the two-photon collisions that we dis- 
cuss do not obey it, producing only dileptons but not 
photons, and the question is how important are those in 
the kinematical range at hand. 



II. THE FORMALISM 

We use the Equivalent Photon Approximation 
(EPA) [T2I [13] to determine the differential cross-section 
for the production of dileptons in Au-Au collisions. Ac- 
cording to this method the effects of the electromagnetic 
fields from the moving nuclei can be replaced by the 
equivalent photon spectrum 
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Following [13] we see that the main contributions to the 
cross-section come from the regions where gi and 92 are 
small (of the order of Wi/7). If both of the momenta are 
small then the total transverse momentum would also 
have a small value and we are interested in studying the 
dilepton production for total transverse momentum up to 
about 0.7 GeV. This is why we work in a semi-coherent 
approach, in which from one of the nuclei we will get a 
coherent electric field, which will correspond to a photon 
with small transverse momentum, while the momentum 
from the other photon can have greater values. This 
means that in this case we won't be getting a coherent 
field from all the nucleus, but that the protons that 
compose it can have an individual effect. For this case, 
instead of using the form factor for a continuous charge 
distribution we will use the one coming from considering 
that the nucleus is composed of Z point particles. 
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where Zi is the number of protons in the nucleus, F(g^) 
is the form factor of the nucleus charge, g^j, is the 
transverse momentum of the photon and Wi is its energy. 
The differential cross-section for the gold-gold collision is 
then given by the product between the photon spectrum 
of each nucleus and the cross-section for the production 
of dileptons from a 2-photon collision: 
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which can be written in terms of the total transverse 
momentum Q = ql -\- q2 and integrated over q2 to give: 




FIG. 2: (Color online) Dilepton production from a semi- 
coherent process. 



For the case when q2± <^ Q± we can approximate 
Q± ~ <l2± ~ Q± so that we can take all the terms with 
Q± from the integral, to get: 
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The factor 2 in front comes from summing the two cases: 
when gi_L is small and when q2±_ is small. Now, using 
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with two parameters, the nuclear radius R (6.55 fm for 
Au)and the width of the nuclear edge which is typically 
about a — 0.5 fm 16^. Starting from this charge distri- 
bution it is not possible to get an analytical expression 
for the form factor, but the integrals of the fourier trans- 
formation can be done numerically, to get a form factor 
of the shape seen in FIG. |3] 
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where Wi = \f{Q\_ + qlz), ^2 = |g2z|and 
rrit = a/M^ -f Q\, we make a change of variables 
from the photon longitudinal momenta qi^ and q2z to 
the invariant mass M and the rapidity y. Then, putting 
2/ = we get: 
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Finally taking the integral over the invariant mass M , we 
get the cross-section as a function of the total transverse 
momentum Qj. and the rapidity y. 
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Since we work in the semi-coherent approach the cross- 
section a^^ is calculated using qi = (wi, Q±, qz) and 172 — 
{qz,0,—qz)- This gives as a result aj^{M,Q±,(l)i,0i). 
The angle 0q is integrated over 27r. The PHENIX de- 
tector covers \ri\ < 0.35 and a total of 180° in azimuth, 
but the data has been acceptance corrected to include 
electrons and positrons from all directions. The only re- 
striction that we must impose is due to the single track 
acceptance condition that p± > 0.2 Ge\^. 
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FIG. 3: (Color online)The square of the form factors plotted 
on a logarithmic scale. The (blue) dashed line corresponds 
to smooth Woods-Saxon charge distribution, the (red) con- 
tinuous line corresponds to resolved discrete protons (but not 
quarks), as explained in the text.) 

As it has been stated previously, we are working in a 
semi-coherent approach. This means that while one of 
the photons is soft and thus sees the nucleus as a uniform 
charge distribution (with the Woods-Saxon shape), the 
other can have a large transverse momenta and thus 
resolve individual protons. For this later case we will use 
the picture of instantaneously frozen nucleons, which 
just means that at any given moment the protons are 
randomly distributed in the nucleus according to some 
weight and frozen in these positions Xm, where x is in 
the direction where the momentum Q± is directed. So 
the form factor can be written as: 
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The charge distribution of the nucleus can be well 
parameterized by the Woods-Saxon expression 



In the amplitude we have the square of the form factor, 
so what we need is: 
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In Figplwe plot this quantity and compare it to the form 
factor from the Woods-Saxon distribution. We see that 
for small momenta they are almost the same but as the 
momentum is increased the Woods-Saxon one goes to 
zero while the one for the point charges doesn't. This 
happens because at small momentum the nucleus is seen 
as a whole and all the protons are acting coherently. 
From eqn 11 we see that for k{xm — a;„) <C 1 the square 
of the form factor |_F(fc)p ~ 1, which means that the nu- 
cleus acts coherently for momenta up to about k ~ 1/R 
but as the momentum increases, the contributions from 
the different protons start to add up incoherently until, at 
high momentum, all the different phases cancel out and 
we are left only with protons interacting individually with 
other protons. So the form factor for the Woods-Saxon 
charge distribution, or for any other uniform distribu- 
tion which may represent the nucleus, is only valid for 
low momenta but stops working as the momentum in- 
creases, due to the fact that the nuclear substructure can 
be discerned at high enough momenta. 



IV. RESULTS 



for small pair transverse momentum Q±, for invariant 
masses smaller than OAGeV . 

In our calculation, if the electron mass is neglected, the 
77 — >■ e~^e~ cross section has coUinear singularity, when 
the leptons move along the directions of the photons. The 
experimental condition ( 13 1 supersedes the cutoff coming 



from finite electron mass: thus we apply the same single 
track acceptance condition in our calculations. 

Since we work in a semi-coherent approach, we keep 
the transverse momentum of one of the photons small, 
with an upper bound given by min{Q±, 1/-R), where the 
whole nucleus acts coherently. In principle, the lower 
bound on q2± is extremely small, related to the electron 
mass, and ultraperipheral collisions with very large im- 
pact parameters would be included in the result. How- 
ever, as we mentioned earlier, these are not included in 
the PHENIX data, so in order to compare it with our 
results we had to restrict the minimum value of q2i_ ■ For- 
tunately, the integral over this variable is logarithmic, so 
we need only order-of-magnitude estimate of its lowest 
value. Roughly speaking, PHENIX only includes colli- 
sions in which two nuclei directly touch each other, or 
b < 2R. The strongest and the weakest electric fields 
of one nuclei inside the disc of the other corresponds to 
distances from R to to 3R. The latter distance we thus 
associate with the cutoff on the minimal q2± value. This 
is needed only for the "softest" photon: the other one 
has large qi± , by kinematics, and thus corresponds to 
field fluctuations inside the disk of the nucleus 1. 



In order to present our results in a way that allows 
them to be compared with data from RHIC, we will 
study the invariant yield which is given by: 
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where atotai ^ 47ri?2 ^ 1.4 . IQ^ [GeV-'^]. The data from 
PHENIX correspond to a minimum bias situation, so the 
collisions considered have centrality in the range — 92% 
which corresponds to Npart = 109 [T7]. Note that very 
peripheral, elastic and diffractive events, both in terms 
of pion production and nuclear physics, are not included: 
however those still can generate dileptons. 

It is also important to notice that the "acceptance cor- 
rected" simply means the geometric acceptance of the 
pair. Separately from this, the detector has a single 
track acceptance condition for the transverse momentum 
of each lepton 



P± > Prnin = 0.2 GeV 



(13) 



needed for them to reach the detector in the current 
magnetic field. The reported cross section [TO] still is 
under this condition. This translates into a truncation 
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FIG. 4: (Color online)Yield versus total transverse momen- 
tum. From top to bottom: coherent, semi-coherent, totally 
incoherent. 

In FigH] we show our results for the yield as a func- 
tion of the total transverse momentum Q±. As expected, 
we see that for small momenta, while it is in the co- 
herent regime, the differential cross-section is large, but 
when the momentum is increased the yield decreases. We 
compare the semi-coherent case with a totally coherent 
case, in which we took the nuclei to be point particles 
of charge Ze, and with the totally incoherent case, when 



we ignore all interferences and consider collisions of pro- 
tons each on each. We see that the semi-coherent case 
lies in between. It starts, for low transverse momentum, 
overlapping with the totally coherent curve and as the 
momentum increases it drops, but it never reaches the 
incoherent curve because we let one of the photon trans- 
verse momentum be small, so that one of the nuclei is 
always giving a coherent contribution. 
For different mass bins the shape of the yield as a func- 
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FIG. 5: (Color online)Yield versus total transverse momen- 
tum for different invariant mass ranges. Open (red) symbols 
are the two-photon contribution, compared with PHENIX 
data (black closed points, taken from Figfl|. From top to bot- 
tom: M=300-500, 500-750, 810-990 MeV. For the lower mass 
range the single track acceptance was relaxed from p± > 0.2 
to p± > 0.1 in our calculations. 

tion of the transverse momentum is the same. In FigjS] 
the yield for three mass ranges is presented and it can 
be seen that the lowest the invariant mass the higher 
the yield, as expected. In order to include the upper 
plot, which is in the mass range 300 — 500MeV, the 
single track acceptance constraint has been relaxed to 
p± > O.lGeV for this case. We compare our results with 
experimental data from PHENIX, which is given by the 
filled black dots in the plot. We see that the contribution 



from the semi-coherent production of dileptons is about 
two orders of magnitude smaller than the experimental 
results. From this we conclude that in the momentum 
ranges explored the dilepton production from two pho- 
ton collisions doesn't contribute significantly to the total 
production of electron-positron pairs in Au-Au collisions, 
but this mechanism is important for smaller transverse 
momenta, when the two photons are a result of the co- 
herent interaction of the nuclear electromagnetic fields. 

The applicability of the EPA requires that the photon 
transverse momenta are small in respect to all other in- 
variants. Therefore, in FiglS^we stop our curves when Q^ 
reaches the dilepton mass squared, where this condition 
is not fulfilled. Note however that in the region where 
the two-photon production has chances to be observed, 
this condition is rather accurate. 

For completeness, we also calculated a contribution of 
this semi-coherent approach for the LHC energy range, 
where 7 ~ 3400. In the unrestricted case, when the 
photon with small transverse momentum is in the range 
< q2±<miniQ^.i/R), the integral over q2± gives the 
usual /n ^''^^'^"'"^ ), so the increase in gamma implies a fur- 
ther "ultraperipheral" enhancement of the process. Us- 
ing the ALICE detector acceptance of |?7| < 0.9 with 
full range in azimuth and a single track acceptance of 
p± > 0.1 [IH] and allowing q2± to be integrated in the 
region just described, there is an enhancement of one or- 
der of magnitude in comparison to our results for the 
PHENIX acceptance and with q2± in the same range. 
However, when we restrict q2± to be between l/3i? and 
min{Q±,l/R) and in this way don't consider ultrape- 
ripheral collisions, the results for the yields that we cal- 
culate for PHENIX and ALICE are very similar and the 
small difference between them (about a factor 2 for small 
transverse momentum) is due to the greater acceptance 
of the latter and also to the different elements used (Z=79 
for RHIC and Z=82 for LHC). 



V. RESOLVING QUARKS 

To end this study of dilepton production in heavy ion 
collisions, we consider the structure of the nucleons, this 
is we resolve partons such that, as before, from one of the 
nuclei we get a coherent contribution (small momenta) 
while from the other one we get the effect from partons 
acting individually. To determine the parton contribu- 
tion we must turn to use the parton distribution func- 
tions (PDF's). 

The nucleons are composed of quarks and gluons and 
the probability that a given nucleon contains a con- 
stituent particle with x momentum fraction of the total 
momentum of the nucleon corresponds to Ji[x)dx, where 
the functions fi are the PDF's for the i type constituent 
{i — u, u, d, d, s, s, c, c, b, b, t, t, gluon) for a proton. 

In our calculations we use data from the CTEQ col- 
laboration [TOj and we only consider the three lightest 
quarks: u and d valence quarks and u,d and s sea quarks. 
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TABLE I: (Color online)Nuinber of quarks calculated by in- 
tegrating the PDF's from CTEQ from different Xmin- 



To determine the number of each kind of quark present 
in a nucleon we integrate the corresponding PDF's from 
Xmin to 1, some examples can be seen in table [11 

In order to consider all the partons that are capable 
of emitting a photon of energy w, we take Xmm to corre- 
spond to that energy. Since in Au-Au collisions at RHIC 
the typical center of mass energy per nucleon is of 100 
GeV, a;mm=0.01 implies that all the partons that can 
emit a photon oi w — IGeV are taken into account. We 
consider partons from protons and neutrons, the PDF's 
for both types of nucleon are related by isospin symme- 
try so for the neutrons the up and down quarks are in- 
terchanged with respect to protons. 

We use Xmin = 0.0025 to calculate the number of each 
kind of quark present in a gold nucleus and then we pro- 
ceed as before, working in the instantaneously frozen pic- 
ture. We randomly select the positions of the quarks in 
the nucleus to determine 
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then this quantity is squared and averaged over different 
quark configurations. The result is plotted on the left side 
in FiglHland compared to the previous case when we re- 
solved only up to nucleons. For small momenta the form 



FIG. 6: (Color online)LEFT:Coniparison between form 
factors. Continuous line:resolving up to protons, dashed 
lineiresolving up to partons. RIGHTiYield versus total trans- 
verse momentum for an invariant mas range of 500-750 MeV. 
Upper curve: resolving up to partons. Lower curve: resolving 
up to protons. 

factor is as before, and the nucleus gives a coherent con- 
tribution, but for large momenta, when the value comes 
from the incoherent contributions of individual particles, 
the result obtained from resolving up to partons is about 
one order of magnitude larger than when resolving only 
up to protons, because the number of participating par- 
ticles is increased since we not only consider the valence 
quarks from the proton but also the sea quarks and all 
the quarks from the neutrons. 



In FigJ6]we also present the comparison for the yield 
as a function of total transverse momentum between the 
two cases that we have studied: resolving up to protons 
and resolving up to partons. We see that for momenta 
larger than 0.1 GeV, the contribution to the dilepton 
production process is larger when we take into account 
the nucleon structure, but it is still orders of magnitude 
below the experimental data. 



VI. SUMMARY 

We studied the production of dileptons in heavy ion 
collisions using the two-photon mechanism, in the dou- 
ble effective photon approximation. It is well known since 
1930's [TT] that for low momenta the contribution from 
both nuclei is coherent and in this case the cross section 
for the process is proportional to (Za)"*. Our main inter- 
est was to look at the regime in which this coherence is 
lost. 

Although our original motivation was to explain the 
RHIC puzzles, we learned (relatively early in the cal- 
culation) that the two-photon mechanism unfortunately 
cannot explain any of them. Yet we persisted and com- 
pleted this calculation, for two reasons. 

One is that with relatively minor modification of the 
PHENIX experimental condition - in particular with re- 
duction of the magnetic field - the dilepton from the two- 
photon mechanism would become detectable and maybe 
even dominant. Thus it would be beneficial to identify 
and study those, in next RHIC runs. It would not require 
significant expense, as we speak about quite large cross 
sections. 

Another reason is a theoretical curiosity: what exactly 
happens with the two-photon cross section when the co- 
herence is lost. We found that when it happens, the 
cross section is dominated by a "semi-coherent" regime, 
in which the momentum of one of the photons remains 
small enough to be represented as a coherent field of one 
of the nuclei, while allowing the other photon to have 
larger momentum and resolve individual particles, pro- 
tons or even quarks. We have shown in this paper that 
the semi-coherent approach gives a greater contribution 
than the totally incoherent case, while still allowing the 
dilepton pair to have relatively large total transverse mo- 
mentum. 

We studied two cases, first resolving only up to protons 
an then resolving them into charged partons (quarks). 
By comparing both results we see that the larger the 
number of particles that are resolved, the greater the 
contribution. Unfortunately, such increase- by almost 
one order of magnitude - is not enough to explain the 
total dilepton spectrum observed at RHIC. 

However, there can be similar semi-coherent regime for 
other processes in QCD, especially at Q ~ IGeV. Tra- 
ditionally, few-GeV momenta transfer is seen as a do- 
main of parton model, with its treated via completely 
incoherent PDF's and cross sections for hard processes. 



And yet, even for gluons, people introduced "color glass" 
fields collectively generated by many hard partons [20] . 
The so called "ridge" phenomenon, recently discovered 
by STAR collaboration [^T , is presumably due to early 
local large-scale fluctuation of color field [22,. A colli- 
sion of such fields with the usual incoherent partons can 
presumably generate a QCD analog of "semi-coherent" 
processes we studied above. 
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